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The electrical conductivity of a pion gas 

D. Fernandez-Frail^ and A. Gomez Nicol^ll 
Departamentos de Fisica Teorica I y II. Univ. Complutense. 28040 Madrid. SPAIN. 

The electrical conductivity of a pion gas at low temperatures is studied in the framework of 
Linear Response and Chiral Perturbation Theory. The standard ChPT power counting has to be 
modified to include pion propagator lines with a nonzero thermal width in order to properly account 
for collision effects typical of Kinetic Theory. With this modification, we discuss the relevant chiral 
power counting to be used in the calculation of transport coefficients. The leading order contribution 
is found and we show that the dominant higher order ladder diagrams can be treated as perturbative 
corrections at low temperatures. We find that the DC conductivity (j{T) is a decreasing function 
of T, behaving for very low T as cr(T) ~ jT consistently with nonrelativistic Kinetic 

Theory. When unitarization effects are included, cr(r) increases slowly as T approaches the chiral 
phase transition. We compare with related works and discuss some physical consequences, especially 
in the context of the low-energy hadronic photon spectrum in Relativistic Heavy Ion Collisions. 

PACS numbers: ll.lO.Wx, 12.39.Fe, 25.75.-q. 


I. INTRODUCTION 

Transport coefficients provide relevant information about the nonequilibrium response of a system to an external 
perturbation. In particular, the electrical conductivity measures the response to electromagnetic fields, the DC 
conductivity corresponding to constant fields, i.e., to the limit of vanishing external frequency and momentum. In 
this work we will be interested in the electrical conductivity of a pion gas at low temperatures, as that formed after 
a Relativistic Heavy Ion Collision has hadronized and cooled down below the chiral phase transition. In this regime, 
the medium constituents are predominantly light mesons, whose very low-energy dynamics is described by Chiral 
Perturbation Theory (ChPT) [3,[3I- We remark that ChPT is a low-energy expansion performed in terms of 
with p a typical meson energy or temperature and A^ ~ 47r/,r — 1.2 GeV and it has been successfully applied to light 
meson dynamics and also to describe several aspects of the low-T pion gas iiiiQ 

In principle, ChPT should provide a meaningful expansion for any thermal observable. However, as we will show 
below, transport coefficients like the electrical conductivity are intrinsically nonperturbative and require a suitable 
modification of the naive chiral power counting scheme to account for collision effects in the plasma, as dictated by 
Kinetic Theory. Similar features have been found in the analysis of trans por t coefficients at high temperatures in 
weakly coupled theories [1,0,^3- fact, from elementary Kinetic Theory |lll| we expect the electrical conductivity 
to behave like u ~ jm for a gas of charged particles of mass m with density A’ch and mean free time r ~ l/F 

with F the particle width. Since in ChPT the pion width is formally F = 0(p®) 0 Q this indicates the need for 
considering non-leading corrections. In addition, the presence of nonperturbative insertions, typically of ©(l/F), 
could spoil the naive power counting of perturbation theory, requiring the sum of an infinite set of diagrams to get the 
correct proportionality constant for transport coefficients. This is indeed the case for high T in scalar 0| and gauge 
|l0| theories, where it has been shown that the sum of certain ladder diagrams provide the dominant contribution 
and coincide with the leading order obtained using directly a Kinetic Theory approach (3 . In Lattice Gauge Theory, 
there are also estimates of the electrical conductivity for lower temperatures but still above the QCD phase transition 
M As for other transport coefficients for the meson gas, the shear viscosity has been calculated in miniiiiii 
in the Kinetic Theory framework. 

One of the main motivations of our work is then to analyze to what extent the previous observations apply at low 
temperatures in a pion gas, within the ChPT formalism. The fact that we are not dealing with a coupling constant 
expansion but with a low temperature and energy one will introduce important differences in the analysis of ladder 
diagrams. As a further interesting point, the DC electrical conductivity is directly related to the small energy limit of 
the emission rate of photons from the plasma formed after a Relativistic Heavy Ion Collision da The photon emission 
rate is a very relevant observable, since the photons leave the plasma almost without interaction^ so that they c arry 

important information about the Quark-Gluon-Plasma [T^, Il3 113 113 the hadronic phase [13 113 113 113 I13l ■ 
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Our analysis is in principle applicable to get an estimate of the very low energy photon rate steaming from the pion 
gas. In this context, we will compare with related results in the literature and with experimental information. 

The plan of the paper is the following: in section E] we will sketch the basic formalism used throughout this work, 
in section nni the leading order for the DC conductivity for low temperatures will be obtained, emphasizing the role 
played by the thermal width and the low-T diagram power counting. Higher order diagrams which could be relevant 
will be analyzed in section ITVl where we will show in detail that for low temperatures one can treat those diagrams 
as perturbative corrections in the ChPT framework. The role of unitarity will be explored in section 0 where we will 
also comment on possible consequences of our results in the context of the photon emissivity of the hadronic gas. Our 
conclusions are summarized in section EH We have also included an Appendix with some useful results for thermal 
loop correlators needed in the main text. 


II. BASIC FORMALISM 

We consider a pion gas at temperature T fully thermalized, i.e, with vanishing pion chemical potentials. The 
response of this system to an external classical electromagnetic field can be described within the formalism of 
Linear Response Theory, which gives the total thermal averaged electromagnetic current to first order in the external 
field (Kubo’s formula) as: 


I dVA):,(x')A« (x-x') (1) 

where = d{t) < [J^^(x), J^^(O)] > is the retarded current-current correlator. 

Here we will always treat photons as an external source, i.e, not thermalized with the pion gas. This is consistent 
with the idea that photons (and leptons) escape the plasma almost without interacting, since their interactions are 
weak. This assumption is commonly used when describing for instance dilepton or photon production, where the 
electromagnetic part of the rate is always decoupled from the strong contribution Therefore, and consistently 
with the above, we will retain only the leading order contributions in the electric charge. In other words, we are 
ignoring rescattering of photons throughout the plasma. 

With the gauge choice of the external field Aj); = 0, the spatial current in Fourier space is proportional to the 
classical electric field: 


< >= 




UJ 




where akj defines the conductivity tensor. In position space: 
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We will be interested here in the DC conductivity, defined as the linear response to a constant electric field, i.e., 
independent of time and space. In that case, one defines the real scalar conductivity a as Keakj = gkjO' where: 


= 1 li„. lim 

6 u;—»-0+ |(f|^0+ UJ 


( 4 ) 


where p = 2Im(iA^) is the current-current spectral function (see Appendix). In 0) it has been also used that, 
by Euclidean covariance, {A^)ij (u;,0) = ^(A^)^( a;,0). We remark that the order in which the zero energy and 
momentum limits are taken in 0 are meant to ensure a physically meaningful answer and correspond to an electric 
field slightly nonstatic but constant in space. For instance, in the reverse order static charges would rearrange to give 
a vanishing electric current However, the small energy and momentum limits are sometimes ill-defined as showed 
below and then one has to be careful with their physical interpretation. 

The DC conductivity is given in terms of a retarded correlator. We will work in the Imaginary Time Formalism 
(ITF) of Thermal Field Theory at finite temperature T and calculate the corresponding time-ordered products, which 
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give the retarded Green functions when the external frequencies iun = 27rnT (n G Z) are continued analytically to 
continuous values as iujn —> w + ie |^. We remark that the ITF has the advantage of dealing with the same 
fields, vertices and diagrams as the T = 0 theory, with properly modified Feynman rules. However, one has to make 
sure that the IT expressions obtained after summing over internal Matsubara frequencies are analytic in the external 
frequencies off the real axis. That will be the case for all the quantities considered here, as showed below. Our ITF 
notation and methods for the calculation of transport coefficients will follow closely those in ^3- alternative 
procedure is to work within the Real Time Formalism (RTF) for which the time contours include the real axis |29| . 
On the one hand, the number of degrees of freedom in the RTF doubles, corresponding to the values of the fields above 
and below the real axis and this complicates the diagrammatic evaluation. On the other hand, the results in the RTF 
are already obtained for real energies, so that no analytic continuation is needed and retarded Green functions can 
be obtained from the time-ordered products by simple diagrammatic rules [30) . An example of a transport coefficient 
evaluated in the RTF can be found in for the shear viscosity in the theory. 



FIG. 1: Leading order diagram for the calculation of the DC electrical conductivity in ChPT without collision effects. 


III. THE DC CONDUCTIVITY IN CHPT TO LEADING ORDER 
A. The ideal gas contribution 

Following the strict ChPT power counting, we should then start by considering the lowest order lagrangian 0{p^) 
(the non-linear sigma model) coupled to an external photon. The lagrangian is given in In the imaginary-time 
formalism we just take the same interaction vertices and electromagnetic current with the usual replacements for the 
Feynman rules, i.e, replacing all zeroth momentum components by discrete frequencies —> iojn = 2'KinT and the 
loop integrals by Matsubara sums, i.e., /^ ^ 

The dominant diagram to the current-current spectral function is that in Figure ^ which gives the following 
imaginary-time contribution to the spatial current-current correlation function : 


0 iiv m = e'^T V f _ {2p + q)j{2p + q)k _ 

{2nr[u;l + W + ml][{u;^ + ,.^)^ + \p + ^^ + ml] 

with niTr — 140 MeV the pion mass. For small external momenta, we can drop the q dependence in the numerator, 
so that we are led to the loop integral (I A. 611 with fc = 1, analyzed in the Appendix. With the results showed there, 
one has to this order: 



Pfc(w.M) 


(^^0+,|gl 


o+)= - -e{\q]^ - CO^) 

TT 


Fi [u;y\e-,T] 

l91 


with Fi defined in (IA.14I) . We observe that the small energy and momentum limits give an ill-defined answer, due 
to the product of the vanishing denominator and the step function (see our comments in the Appendix). In other 
words, there are regions arbitrarily close to the origin in the {uj, |^) plane where we get different limits for the spectral 
function: it vanishes for the timelike region < uj'^ and grows with l/\<f\ in the spacelike one \q\'^ > . Therefore, 

in this case instead of specifying a particular order in taking the limits, we consider the Fourier transform of the 
contribution ® back to position space, according to 0 . We take electric fields slowly varying both in time and 
space, so that only the small w and limit of the spectral function is relevant. When doing so, we get a divergent 

contribution for the conductivity a = — dr to leading order in the derivative expansion of the electric 

field, where the dependence with the spatial length of the system is explicitly shown. This is physically expected 
since, as commented in the introduction, we are considering pions with zero width or infinite mean free path. In fact. 
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FIG. 2: 


(a) A two-loop diagram contributing to the pion width. 


(b) Leading order current correlator with dressed pion lines. 


this is nothing but the ideal gas result (no collisions) since by considering only the diagram in Figure ^ we are not 
taking into account the effect of pion interactions. Actually, to this order one would get the same answer, say, in 
scalar QED. Therefore in an ideal gas, the charged free pions travel throughout all the gas without interacting giving 
an infinite conductivity. Thus, for an external field constant over an infinite region, we should consider the effect of 
pion interactions so that the DC conductivity is actually proportional to the (finite) pion mean free path. This will 
require a more careful evaluation of ChPT diagrams, as we will show below. 


B. Beyond the ideal gas: finite pion width. 

Our previous discussion makes it clear that in order to get finite and physically sensible results for the DC con¬ 
ductivity, one has to consider the effect of a finite pion width. In the medium, the pion dispersion relation is 
= rn^ + g{po^ \p\'iT) with g a complex function. The pion width is Tpdpl) = —lnig{Ep, |p|)/(2i?p) << Ep in ChPT 
with Ep = v1?F+ and we will neglect the ChPT corrections to Reg, since they are not relevant for our purposes. 
The lowest order contribution to Img in ChPT is 0{p^) and can be calculated either by Kinetic Theory arguments 
0 or by evaluating the relevant diagrams 0 like the two-loop contribution shown in Figure The detailed results 
for Fp can be found in those references. For most of our purposes here it will be enough to consider the dilute gas 
regime where ns(£'p) « exp(—Ep/T) << 1: 


j^DG _ 
p 


1 

2 


(2^)3'^toRsl 2F;feF;p 


exp{-Ek/T) 


( 6 ) 


where the Mandelstam variable s = {Ek + Ep)^ — \k + pl"^ > 4mJ and crfj) is the total elastic pion-pion cross section, 
that can be expressed as: 


QO_ ^ ^ 

^ E E(2-^ + (7) 

J =0 1=0 

in terms of partial waves of definite isospin / and angular momentum J, for which we follow the same conventions for 
pion scattering as e.g. in 0 To lowest order in ChPT, i.e, |tp == 


TT-TT/ \ _ 


1 - 


16m^ 37mf 
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( 8 ) 


where /^r — 93 MeV is the pion decay constant. Since the integrand in m is weighted by the Bose-Einstein factor, 
we expect that taking the leading order for crfJJ is enough for relatively low temperatures. We will comment below 
on the effect of considering higher order terms, for instance those respecting unitarity. 

The result m shows clearly how the pion traveling along the (dilute) medium with momentum p acquires a finite 
width by collisions with the medium constituents, the factor s{s — Am^)/{2EkEp) bei^ proportional to the relative 
velocity between the particles with momentum p and k, expressed in relativistic form 0. 

We consider then a pion propagator whose spectral function (see Appendix) is: 


Ep ^ {uj 


1 

EpY + Tl 


1 

[u: + EpY+Tl_ 


p^{uj, 1^) = 2ImiG'^(w, 1^) 


(9) 
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Thus, the complex propagator is: 


G{z, \p\) 


1 

E^-[z + zrpSgn(Imz)]^ 


( 10 ) 


for Im z ^ 0, which has a cut on the real axis and is analytic everywhere else in the complex 2 plane. 

We will denote by a double line, the pion propagator “dressed” with a nonzero width Tp, so that the leading order 
contribution to the current-current correlator is now given by diagram (b) in Figure|21 Its imaginary-time contribution 
is now 


n— — oo ^ ' 

where Gjsliujn] = G[iu}n\ in ifTIIIl for n 7 ^ 0 and Gp{n = 0 ) = l/(rp -|- E^) and where the |(^ ^ O'*' limit is now free 
of the ambiguities of the ideal gas case. We will evaluate the Matsubara sum in HU following the same procedure 
as in M- In this case, the result can be expressed in terms of integrals along the two cuts of G at Imz = 0 and 
Imz = —I'm of the corresponding discontinuities of the G function. The result is an analytic function for ii/m off the 
real axis, so that we find for the retarded correlator (A^(a;) = ^ w -f ze)): 


[A«]J(o.,|g1 = 0) = de^y ^|pf {G^(o) [G«(tu) - G^(-o.)] 

+ 1"° + G^iy)] [G^{y + iz)- G^{y - u;)] } (11) 

where we have omitted the |p| dependence in the propagators. Now, it is not difficult to see that the leading order 
term in the previous expression in the limits w —> O'*" and T << Ep is an G(l/r) contribution coming from the 
products G^G^ which in this limit behave like a S function: 

/ OO 

dyG^(y)G^(y)F(y) iF(Ep) + F(-Ep)] (12) 


for Tp << Ep. 

This type of “pinching poles” contribution is crucial for the analysis of the conductivity. On the one hand, it is 
expected from Kinetic Theory considerations, as discussed in the introduction. On the other hand, it will lead to 
consider the effect of higher order diagrams, as we will do in the following sections. Up to this point, our result are 
analogous to those in 11113 before specifying the explicit form of Tp and discussiM the power counting. Similar 
“pinching poles” contributions are obtained in real-time formalism calculations [11 [^ . 

Thus, we get for the DC conductivity to leading order: 

where we have used that ri'g{x) = —nsix) [1 -I- nsix)] /T and rzs)—x) [1 -I- nsi—x)] = 7 Zb(xU 1 -I- nsix)]. In Figurej^l 
we have plotted the leading order conductivity taking for Fp the exact expression given in jj, H including all factors 
of rzB, comparing it with the different approximations considered below. 

The result (tni is very important for our purposes. It gives a finite and well-defined answer for the conductivity, 
with the correct dependence on the pion width and the density as expected from Kinetic Theory arguments. However, 
as commented before, the nonperturbative nature of the DC conductivity, reflected in the l/F dependence, will make 
it necessary to examine higher order diagrams. 

Let us evaluate the previous result for the conductivity in the dilute gas regime, which is valid for the range of 
temperatures of interest here. Thus, we replace the results m and © in m, where we neglect the rz^ term. The result 
is showed in Figure El (solid line) where it is clear that the dilute gas regime in this case is an excellent approximation 
for temperatures below the pion mass. The conductivity of the pion gas grows for decreasing T, approaching the 
ideal gas limit discussed in section nun Note that both the pion density and the pion width vanish exponentially 
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FIG. 3: Leading order for the DC conductivity. The left panel shows the dilute gas mproximation (solid line) compared to 
the exact numerical result obtained by equation llT^ with the full width calculated in P,|3 and to the T « rtin limit (dashed 
line) given by m- In the right panel, the behaviour near T ^ m-K is showed, where the dashed line is given by equation C3- 


for T 0+, so that it was not clear what should b e the r emaining T-dependence of their ratio. As we will show 
later, for very low temperatures T « rriTr, a ^ m-,^ . This behaviour is nothing but the nonrelativistic limit 

predicted by Kinetic Theory arguments and indicates that there could be nonperturbative contributions of higher 
order diagrams in the very low T region, where the typical l/T insertions become larger. 

For this reason, in the next section we will concentrate on the T << ?7i,r limit. As we shall see, in this limit the 
analysis becomes simpler and our power counting arguments will be more clearly illustrated. 


C. Very low temperatures 

In the T « limit, the relevant integrals in the previous expressions are dominated by the very low momentum 
region, weighted by ns(iffe) ~ exp(—TO,r/T) exp(—fc^/2TO,rr). Thus, if we change p ^/m^p in 11311 and k 
y/rriT^Tk in 0, every power of a three-momentum integration variable p/ in those integrals counts as y/rn^ /ttItt = 
yjT/m^ and is therefore suppressed. Taking into account that in this limit s(s — 4m^) ~ 2m7r|p — k\ in 0, the 
angular integration in the expression for the thermal width can be performed and the low-T result can be expressed 
in terms of error functions: 


with 


yLT 
^ p 


TO’ 


Epf^ 


exp{-m^/T)f{p) 


(14) 


f{p) = 


23 

5127r3 


exp{-f/2) + 


/|l^Erf(p/V2)' 

V 2 p 


(15) 


and p = p/y/niTrT. What is important from this expression is that for Ep ^ TOtt, Tp ~ exp(—TO7r/T) times 

a dimensionless function of p only so that we get readily the dependence with T, TO,n.,/7r of the thermal width for 
very low temperatures, as the product T ^ nuafJJ with n ^ (-s/rn^Y (mean density) and v ^ ^/mYTIrriT^ (mean 
velocity) in the nonrelativistic limit of very small pion momentum compared to its mass. We recall that the same 
nonrelativistic behaviour is recovered at very low temperatures for the shear viscosity, which grows with \/T [Hill. 

Replacing the above result for the width in the T « m-^ limit of m yields: 


^( 0 ) 

^LT 




to: 


— / dpj^exp{-p /2) ~ 2.7e — 


(16) 


Therefore, for low T we get a ymY/T behaviour for the conductivity in units of TO,r- This result will be crucial 
for our chiral power counting analysis of higher terms in the DC conductivity and gives the behaviour near the origin 
plotted in Figure l^l where we compare 11611 with the dilute gas approximation. 
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D. Higher temperatures 

For temperatures up to the pion mass we can trust the dilute gas approximation, as shown in Figure |21 On the 
other hand, we do not expect ChPT to work much beyond that range. A fit to the numerical results for temperatures 
around assuming a single power form aT^ gives: 


iT<m^) (17) 

The above behaviour with T can be understood as follows. Consider the high temperature limit T » in our 
previous expressions, which amounts to ignore the pion mass. Thus, for an integrand of the form f{p)nBip), this is 
a better approximation for functions / dominated by their high energy behaviour. In fact, in some cases the high-T 
curve is approached for not very high temperatures, even when the dilute gas approximation still remains valid and it 
is not a bad approximation to neglect the pion mass in the dilute gas expressions, in order to get, at least qualitatively 
the behaviour for T<m,n.. Consider for instance the total density of charged pions N^h = (I/t^"^) / dkk^nsiEk). In 
the T >> ttItt limit we replace nsiEk) — l/(exp(fc/T) — 1) and obtain Nch — ‘2C,{3 )T^/tt"^ = 2AT‘^ whereas 
taking both the dilute gas and high-T approximations, i.e. nB{Ek) — exp(—fc/T), gives Nch — 2T^/7r^, which 
around T ~ TOtt is only about a 10% away from the exact curve. Something similar happens with the pion width 
Fp, for which the high-T limit in the dilute gas expression m with the lowest order cross section in (ISJ gives 
Fp ~ 5T"‘p/(127r^/^). Note that it is crucial for this argument that the cross section grows with energy which in fact 
is not the physical case, since unitarity bounds the partial waves. Therefore, we see that implementing correctly the 
unitarity behaviour may change the behaviour of the thermal width and hence of the conductivity for intermediate and 
higher temperatures. We will examine this aspect in section[3 Ignoring those possible corrections, the mean thermal 
width f = J TpHBiEp)/ f nB{Ep) ~ 5T®/(47r^/))) in this limit, in agreement with the analysis in 0. Therefore, 
we expect qualitatively for the conductivity a ~ e^A"c?i/(w 7 rf') ~ e^m 7 r(/ 7 r/w,n.)‘*(m 7 r/T)^ as our above result 11711 
roughly confirms. Performing the same approximations, i.e., dilute gas in the chiral limit, directly in the expression 
(tT^ gives = 27re^/^/(5T^) which is worse than the fit (tTTIl at T<m 7 r, mostly because the integrand in 0 grows 
much slower with p than in the case of the width. 

We also note here that at higher temperatures, i.e, closer to the chiral phase transition, one has to be careful 
with the approach we are following throughout this work, since the thermal width becomes comparable to the mass, 
as a consequence of the loss of validity of ChPT. In fact, using our previous high-T expression for the width gives 
f ~ m^/25 at T = 100 MeV but f ~ 1.2m^ at T =200 MeV. 


IV. HIGHER ORDER DIAGRAMS 

Our previous analysis shows that the leading contributions to the DC electrical conductivity come from nonpertur- 
bative 1/F insertions for every pair of pion lines carrying equal momentum as the external four-momentum vanishes 
(double lines). Therefore, it is not obvious that diagrams of higher order in the chiral power counting are still 
suppressed with respect to the leading order analyzed in the previous section. 

We recall that according to Weinberg’s power counting, the chiral power of a given diagram with two external photon 
lines is e^p^ with D = 2Nl + J2d=2i^ ~ Nl being the number of loops and Nd the number of vertices coming 

from the lagrangian of order p‘^ in derivatives and meson masses. However, the presence of “pinching poles” modifies 
this picture. Let us assign a “nonperturbative” factor Y to the contribution of every two internal double lines to the 
conductivity, so that we define = e^triT^Y and we call X = OijP') a chiral power, typically X ~ . 

Thus, the diagram in Figure[5|D analyzed in the previous section is 0{Y) instead of the 0(X) assigned by Weinberg’s 
power counting. 

The next task is then to identify those diagrams which are dominant with this new power counting. In the case of 
a scalar X(j/ theory, one has @ T = 1/A^ and the dominant contributions are ladder diagrams of the type showed in 
Figure 01 which are all 0(1/}?) and have to be resummed, giving the results expected from Kinetic Theory. Other 
diagrams which naively could give larger contributions, like the “bubble” diagrams showed in Figure 0] have been 
shown to be subdominant. The same conclusion about the dominance of ladder diagrams holds in gauge theories 

illSIIl. 

The arguments based on the topology of the diagrams in the scalar theory apply directly to our case, so that we 
are led to examine the same class of diagrams. However, it must be pointed out that in ChPT there are important 
differences: first, we are actually dealing with a double counting in X and Y as explained before, which is a consequence 
of the importance of chiral loops in the power counting. Second, we are allowed to have derivative vertices and multiple 



FIG. 4: Ladder diagrams contributing to the electrical conductivity. 



FIG. 5: A generic bubble diagram with fc + 1 bubbles (fe > 1). 


vertices coining from lagrangians of different orders. And finally, there are two types of particles, charged and neutral, 
running in the internal lines and loops. This motivates our detailed analysis below, which main conclusions will be 
that here the ladder diagrams also dominate over other diagrams like the bubble ones, but their influence is not as 
important as it could be and can be safely considered subdominant with respect to the leading order result, within 
the ChPT validity range. 

Let us start by considering then the ladder diagrams showed in Figure0]where all the vertices come from the d = 2 
lagrangian. Note that only the internal lines running along the ladder can carry equal momentum when the external 
momentum vanishes, so that we keep the double line notation for those lines and the single line for the rung loops. Our 
power counting gives then 0{X^Y^^^) for the diagram with n rungs. Therefore, those diagrams could in principle 
be of the same or even higher order than the zeroth order diagram analyzed in the previous section, depending of the 
relative size between X and Y. Now recall that according to our previous analysis, the nonperturbative Y factors are 
bigger the smaller the temperature, where typically Y ~ ^Jrr^^jT. It is therefore at very low T where we have to be 
more concerned about these contributions. The detailed analysis in this regime is done in the next section, where we 
will see that at low T one has to account also for the Bose-Einstein distributions emerging from the loops, which will 
give rise to additional suppressions. 


A. Ladder diagrams at very low T 


For the analysis of the generic ladder diagram showed in Figure 0] we will follow similar steps as in M- The 
contribution of the diagram with k rungs to the current-current correlator is now: 




(zt'm,0)=e^r ^ 


rW/- , ■ ■ 

,, 2pi • AV(JW„ + IVjn^liOn^Plj 
[ZTT)'^ 

xGf3iiU!n, \pi\)G(jiiuJn + iv-m, |pi|) 


(18) 


where is the effective 7r+7r“7 vertex made of k rungs, as depicted in Figure 01 where the sum is made over all 
possible ways of inserting internal lines and vertices corresponding to charged and neutral pions, respecting charge 
conservation and including the corresponding combinatoric factors. Note that the two pairs of lines attached to 
the external lines are always charged, denoted in the figure by -I-, — symbols. The lowest order effective vertex is 

The effective vertices satisfy a recurrence relation which ultimately leads to integral equations corresponding to 
linearized Boltzmann equations I1II3 The same equations for the effective vertices are obtained in the real-time 
formalism In our present case, the integral equations become a little more involved due to the presence of 

neutral and charged particles, as well as derivative vertices. However, for T « m,r an important simplification takes 
place: The leading order of the sum in Figure ^ for a given number of rungs k is proportional to a single diagram 
where all vertices are replaced by constant vertices //J and therefore the spectral function of the rung integral is 
given by that of the loop correlator Jq in (IA.15II . analyzed in the Appendix. 

Let us justify our previous statement. We have seen in the previous section that the presence of a pair of double 
lines forces its momentum to be on shell at pq = ±Ep. Thus, consider the one-loop subgraph made of two “incoming” 
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double lines and two “outgoing” ones. The two vertices contain different combinations of powers of two momenta or 
two masses in each vertex, where the momenta can be “external” (double lined) pi , P 2 or “internal” k and k + pi ± p 2 
(see details below). In addition, as we will see below, the integrals over pi and p 2 contain Bose-Einstein factors that, 
according to our very low T counting, suppress powers of pi, 2 /™ 7 r- Therefore, the leading contribution for T « 
in or p\ 2 vertices is proportional to m‘^Jo{k), since Ei ^2 — As for the powers involving “internal” momenta, 
using that Ei + E 2 — 2 m.,r >> |pi — ^| >> Ei — E 2 = C>(|pli, 2 /^ 7 r) and the relationship among the different 
one-loop thermal integral with momentum powers in the numerator discussed in Appendix A of 0 , we have checked 
explicitly that we are left, to leading order, only with rn^Jo, m^Ji or m’^J 2 contributions to the imaginary part 
(spectral function) of the loops, which, as our analysis below will show, provides the relevant contributions to the 
ladder diagrams. The imaginary part of the Ja integrals is given in the Appendix. For T << TOtt we can replace the 
Bose-Einstein distributions by Boltzmann exponentials and the relevant integrals can be explicitly evaluated, with the 
result that both for the thermal and for the unitarity cuts, the leading order of Im are proportional to m^Im 
for a = 1,2. 

Therefore, for T << TOtt we have 


^'■ch — 



4fe 


(19) 


with a purely numerical factor and 


AW(i 


lUJm \ 


^ X! ^\iuJn+iVm,iuJn,P2) 

Gi3{iuJn, \p2\)Gf3{iU!n + Wm, |P2 | ) 


Uj' - i{uJn - OJm) 


( 20 ) 


where p(uj, |^) = 2Im J(f (w, |^) = 2Im/^(tu, |^; 0) according to our notation in the Appendix. The above recurrence 
relation is depicted in Figure|^and allows us to proceed along similar steps as in [lO). By induction, we have that the 
only singularities of + z,p) are the same two cuts as the propagator product in (EHll . i.e, at Imx = 0 and at 

Imz = —Vm- This allows to perform the frequency sum in Hi in a similar manner as we did in the previous section 
and also the sum in the equation for the effective vertex CT . where in the latter case we have to consider, in addition 
to the cuts contribution, the pole contribution dX z = ui' + iuJm- As in the previous arguments on the analytic 
structure of the effective vertex imply that the analytic continuation Wm ^ w -fie is well defined. Altogether, we find 
for the low-T conductivity to leading order in l/T, using itT^ : 


with 


dfe) _ 


g 2 Q,(fe) 

ST 



^^^nB{E^)[l + nBiE^)]A('^\\p,\) 


( 21 ) 


A«(|pi|) = ^^^^A(>^-^H\P2\){[nB{E2-E^)-nB{E2)]piE2-E,,\p2-Pi\) 

- [nB{E2 +El)-nB{E2)]p{E2 +Ei,\p2+pi\)} ( 22 ) 

where we have also used that p(w') = —p{—u}') and we have denoted for the effective vertex + ie,x — ie,p) = 

2pA^^'> {\p\-, x) and = {A^^\\^-,Ep) + A^^\\p\-,—Ep))/2 so that = 1. For simplicity, we have denoted 

Fi ^2 = Epj^^p^ and so on for ri _2 and omitted their explicit dependence with \pi\, |p 2 |. 

Let us comment on the result 123 . The two terms in the right hand side correspond to evaluate the imaginary 
part of the rung loop at pi ± p2 for on-shell pi^2 “external” momenta. That is, those contributions arise from the 
s,t channels of the tttt elastic scattering amplitude of the double lines in Figure n to one loop. Here we call by 
convenience S = (Ei + E 2 ,\pi +P 2 I) and T = {Ei — £' 2 , |pi — P 2 |)- This interpretation will be very useful in our 
analysis of the different ChPT diagrams. In fact, note that the remaining u-channel contribution would be given by 
the contribution to the effective vertex of the bubble diagrams in Figure |31 since the loop integral of the bubble does 
not depend on pi and P2, only on the external photon momentum q, i.e, we can identify U = {uj, ()). Now we remark 
(see Appendix) that the imaginary part p{E, jQI) is nonzero both for E^ > jQp + 4to^ (“unitarity” cut) and for 
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FIG. 6: Recurrence relation between ladder vertices. 


< IQP (“thermal” cut). Since s = > 2m^, the s-channel is given by the unitarity contribution in fOl) with 

E = El + E 2 ,Q = Pi + P 2 , while t = < 2m^ and thus the t-channel is given by the thermal part in (IA.12II with 

E = El — E 2 ,Q = Pi — P 2 - Note that a consequence of this analysis is that, as announced, the bubble diagrams in 

Figured are indeed suppressed in the w ^ O'*" limit since u = U^ = uj^ does not fall within any of the two cuts. This 

will be confirmed in section iwn 

At very low T, there are further simplifications of the above equations, namely n'^{Ei 2 ), n{Ei + E 2 ) « n{Ei 2 ) — 
ex'p{—Ei^ 2 )/T. Consistently, in the 5-channel, we can neglect the T-dependent part in the unitarity contribution iTOll . 
In the T channel, we can approximate in iTOHl the logarithm in the right hand side by e (l ~ ^ = 

_ ^ 2 !). Since yo > tUtt, the leading contribution in that term comes from the n_B(i ?2 — Ei) 
term in (P|). giving a net 0{e contribution, the same as in the 5-channel and in the rir 2 in the denominator. 

The latter approximations are valid in the dilute gas regime, which as commented previously, extends to 
In addition, for T « in the 5-channel, the product of the Bose-Einstein distributions n{Ei)n{E 2 ) implies 
El 2 — ViiTr + [|pi 2 p/( 2 TO,r)] and hence tT 2 p — \Pi ~ P 2 |/( 2 TO,r)- Finally, we have for T << rriT^'. 


LT 3 


A 

m-n 


A 


(k) 

LT 


(^\/m^Tyi'^ 


1 


(m. 

~T 

(Pm 


(Pyi 


_ Vl Ak) 

(2^)3 /(yi) 


{^•s/m^Tyi 


exp 


yl 


327ryf J {2tt)^ \yi - y2|/(?/2) 


. At ^\\/'m^Ty 2 ) exp (- ^ 


\yi - ml 


exp 


2 2 1 ivi ~ yir 1 1 ^ ^ i2 

yi E y2 — T tv: -T-|2-T yi — 2/2 

2 \yi - 2 / 2 r 2 


(23) 


(24) 


for fc > 1 , where the function / is defined in OSl and we have rescaled pi ,2 = •\/m)(Tyi ,2 ( 2 / 1,2 = |yi, 2 |). 

From the previous equations CT-(I^ we can draw one of our main conclusions: the low-T effective vertex 
(VwjrFyi) is T-independent and hence, from 12311 the correct order for the contribution of the fc-rung lad¬ 
der to the conductivity in our power counting scheme is = e^mTrO{X^^/rrh^jT) = e^mT^O{X^Y) rather than the 
0{X^Y^'^^) following from direct inspection, as discussed at the beginning of this section. The independence of the 
LT effective vertex on T follows from lOl) by induction, since = 1. Then, the only dependence with T in is the 
1/-\/T, the same as the lowest order CSl. We can also understand why the “direct” counting misses k inverse powers 
of Y. The reason is that the contribution from the spectral functions introduce an extra -^TIm-^ factor inside the 
integrand in |23Ji. These factors come from ps ~ |pi — _P 2 |/w,r and pr ~ El\pi — P 2 I before rescaling the integration 
variables. This reflects the fact that the imaginary parts are small at low T due to the presence of the Bose-Einstein 
factors and represents one of the main differences between our analysis and the standard one in a coupling constant 
perturbation theory. The remaining O(X^) contribution comes from having k chiral loops. To check explicitly our 
conclusions, we have evaluated numerically the one-rung contribution /c = 1 in setting = 1. We obtain 

— 0-06. The reason why this gives a slightly higher value than expected by power counting is the presence of 
a collinear effect in the T-channel, coming from the region yi ~ 2/2 in (EH). The integral is convergent but that region 
enhances its value with respect to the 5 channel. In fact, for k = 1 the T-channel contribution is about five times 
larger than the 5-channel one, which has the expected chiral reduction mentioned above. Collinear effects of this 
kind are characteristic of these analysis and play a crucial role in gauge theories mill. In our case, this numerical 
enhancement does not spoil the perturbative nature of these contributions. 

We can then conclude that the ladder diagrams at very low T, although representing the main contribution beyond 
leading order, are perturbative in the ChPT scheme. In numerical terms, the very low T conductivity is given by itTCI) 
with a theoretical uncertainty of only about 5% in the multiplying coefficient. 
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FIG. 7: Estimation of the contribution of ladder diagrams for higher temperatures, with and without derivative vertices. The 
points correspond to the exact expressions without using the dilute gas approximation. 


B. Ladder diagrams for r<m,r 

Our previous analysis shows that the ladder diagrams are still chirally perturbative in the very low T region, 
even though we should not miss the fact that their correct counting gives a much larger contribution than their 
standard ChPT one, namely = 0{p^^) instead of the 0{p^^) given by Weinberg’s theorem. The reason 

is the presence of the nonperturbative Y factors, as we have seen. Following our chiral counting arguments for 
temperatures of the order of m,r, we realize first that now Y ^ 0(1) around T ^ m,r- This means that we expect 
again for the fc-ladder diagram, coming from the rung loops. It must be reminded though that 

for these temperatures, the simplification of reducing any ladder diagram to one with constant vertices does not hold. 
In fact, for higher temperatures we expect derivative vertices to become increasingly important, but as long as T 
remains within the ChPT applicability range (not much above the pion mass) we can consider all those diagrams 
perturbatively suppressed. 

To check our previous comment we have evaluated numerically the ladder diagram with one rung and constant 
vertices, setting in PT |) - (P|l k = 1 and = 1 and performing only the dilute gas approximations mentioned in 
the previous section but not those specific of T << 77i,r- We have also considered, for comparison, the contribution 
of a ladder diagram with one rung but with one possible combination of derivative vertices, namely, that obtained 
by taking fc = 1 in PT |) -(P |) but replacing the constant factor by |pi|^//^ in the integrand. The results 

are showed in Figure [7| where we have showed for comparison some points calculated without using the dilute gas 
approximation. Clearly, the derivative vertices become the more important ones at moderate and high temperatures, 
showing that the ladder diagrams have to be summed at temperatures close to the phase transition, where the ChPT 
power counting fails. 


C. Other diagrams 

Following our previous arguments, it is not difficult to conclude that, as far as the ladder diagrams analyzed in 
sections IIV Al and IIVBI are perturbatively controlled in the chiral expansion, so are the same diagrams with vertices 
coming from the lagrangians Cd with d > 2 and also ladder diagrams with rung loops made with more than four pion 
vertices. 

As for the bubble diagram showed in Figure^which could in principle give nonperturbative contributions of 0(1"”), 
our previous discussion based on the scattering of internal lines has already suggested that these diagrams do not 
receive contributions of 0{1/T) to the conductivity. Let us show this in more detail. First, we define a /c-bubble 
effective vertex with k>l (see Figure El exactly as we have done in CHI for the ladder vertices A^^^. Given the 
structure of the bubble vertices, one readily realizes that, unlike the ladder ones, + is independent 

of the frequency and momentum of the external pion leg if only constant four-pion vertices are considered. 

Therefore, the contribution of those diagrams to the conductivity vanishes by parity, for instance by making pi —> —pi 
in the integral analogous to HU. Thus, one needs derivative vertices to yield a nonvanishing contribution. In fact, 
following the parity argument, the only contributions of a fc + 1-bubble diagram which do not vanish by parity are 
those where all the four-pion vertices contain the factor pi ■ Pi+i {i = 1 ,..., fc) where pi is the momentum running 
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through the i-th bubble. We remark that we are considering only pion vertices coming from £ 2 , higher order vertices 
being suppressed by Weinberg’s power counting and the same holds for bubbles made with more pion lines, say with 
vertices with more than four pions. Thus, the bubble effective vertices satisfy a recurrence relation of the form: 


J Jtt 


n— — oo 


(25) 


In fact, we should take into account that here can be also neutral pion loops in bubble diagrams with more than 
two bubbles, implying different vertex and combinatoric numerical factors in front of pi ■ Pi+i for the Tr'^Tr'^Tr+Tr”, 
7r'''7r“7r'''7r“ and vertices. Thus, we really have to define two effective vertices 77r'''7r“ and 77r°7r°, the latter 

starting with one bubble and one has then two coupled recurrence relations between those vertices, instead of the single 
one in (1^ . However, that merely complicates the calculation but does not change our main conclusion, namely that 
bubble diagrams are subdominant. Hence, for simplicity, let us show it for the case in which all dvr vertices contribute 
simply as pi ■ pi+i/f^ and all internal lines are charged. Proceeding then as in section ITV Al we realize that the 
calculation is particularly simple now since the vertex (|23 is still independent of the frequency of the external leg 
ium and therefore there is no pole contribution, unlike Then, the Matsubara sum is just the same as that in 

CD- We get, after analytic continuation ii/m uj + ie and to leading order in l/T and w —> O’*": 


p) = - / ^^2 P- p-2 n{E2) [1 + n{E^)] 


(26) 


for /c > 1, where the last line holds by Euclidean covariance and defines the scalar function F, which satisfies then 
the recurrence relation: 




iioKT) 


F^'^-^\uj;T) = 2HT) 


'iujI{T)Y 


(27) 


where I{T) is the integral appearing in the lowest order conductivity in lldll . namely = Ae^I{T) and we have 
used that 73^°^ = A® = 2pj so that = 21 (T). 

From the above results, we note that the two main differences between the bubble effective vertices and the ladder 
ones are, first, that for small w, the contribution of the k + 1-bubble diagram to the conductivity is O{oj'^) and, 
second, that the bubble effective vertices are not all real, is real for k even and purely imaginary for k odd, to 
leading order in l/T. This confirms our discussion in section H V Al about the role of the different diagrams as related 
to the “scattering” of the internal pion lines. The one-bubble diagram amounts to the tree level contribution to 
the scattering in ChPT which is always real. The first diagram giving a nonvanishing contribution to the scattering 
spectral function is the two-bubble one, but it vanishes as oj ^ 0“*" as corresponds to the U channel. 

The simple form of the bubble vertices allows to perform explicitly their contribution to the conductivity, since they 
form basically a geometric series according to m- The fact that bubble diagrams can be summed is a well-established 
result and holds also in other analysis of transport coefficients Q- In fact, we remark that in our case it is consistent 
to sum them, following our power counting arguments, since the bubble diagrams are neither reduced by chiral loops 
like the ladder ones (all lines are double) nor by higher order pion vertices if all pion vertices come from £ 2 - 

Thus, the contribution of the sum of the bubble diagrams to the conductivity is of the form: 


^bub 


lim — 
u;^0+ T 


(Yp n{Fp) [1 -I- n{Fp)] 

(^ 


Re ^2p- 5 ^( 0 ;,:^ 

fc=i 


2e2 lim Re VFW(a;) = a(°)(r) lim Re 

o[P)P°\t) 


(28) 


and the same conclusion, i.e, (Jhub = 0{uY)P^\T) is reached when considering the full contributions of charged and 
neutral diagrams, as commented before. We wish to stress that the leading order in the bubble vertices comes from 
the “pinching poles” in products GrGa and carries the 0{u!/T) contributions analyzed above. Taking the next to 
leading order could give a nonvanishing contribution as w —*■ 0"'' but with no l/P terms (the I(T) integrals above) 
and therefore suppressed by the ordinary chiral power counting of the loops so that, typically, the ujI{T) factors in 
the sum (1^ would be replaced by an 0{pY contribution. Therefore, our conclusion is that the bubble diagrams are 
subdominant for all temperatures within the validity range of ChPT. 
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V. UNITARITY, HIGH T BEHAVIOUR AND PHOTON PRODUCTION 

One of the conclusions of the preceding sections is that the behaviour of the conductivity with temperature is 
greatly influenced by that of the thermal width, as expected from physical arguments. At the same time, the thermal 
width in the dilute gas regime @ grows with T according to the growing of the partial waves and cross-section with 
energies. We have been using the low-energy leading order in ChPT 0 which grows with energy s. However, we 
know that partial waves should be bounded in energy, as a consequence of unitarity. We recall that the unitarity 
condition for the S'-matrix translates into 


= cr2p(s)|t/jp (29) 

for the elastic tttt scattering partial waves (s > and below any inelastic threshold) where a 2 p{s) = -^1 — 4to ^/s 
is the two-pion phase space. The unitarity condition implies that the modulus and the phase 5ij (the phase 
shifts) of the partial waves are related as \tjj\ = sinSu/a 2 p which means in particular that |t/j| is bounded as energy 
increases, as commented above. 

The ChPT partial waves only satisfy unitarity perturbatively, for instance = cr 2 p\t^^]\'^ where the chiral 

expansion is tjj = -|- .... They are not bounded and they cannot reproduce a resonant behaviour, like 

those appearing in the / = J = 0 (the a or /o(600)) and / = J = 1 channels (the p(770)). However, it is possible to 
construct amplitudes satisfying exact unitarity and matching the low-energy ChPT expansion. In doing so, we will 
follow the so called Inverse Amplitude Method (lAM) relying on its remarkable success to describe meson-meson 
scattering data up to 1 GeV and to generate the low lying resonances . We remark that through this unitarization 
procedures, physical resonances are generated dynamically with masses and widths in agreement with those listed in 
the Particle Data Book . Furthermore, the unitarization procedure can be extended to finite temperature, providing 
a consistent description of the thermal properties of resonances, such as the widening of the p in the thermal bath 0. 
The lAM unitarized partial waves for tttt elastic scattering are simply given by The same 

formula is valid at finite temperature, where the fourth order thermal amplitudes have been calculated in [^. 

We are then interested in the effect of unitarity on our previous results. The first important comment is that 
unitarity is not expected to change our low T results, where only the low-energy region of the tttt scattering contributes. 
We expect more significant corrections for T near m^n-. The effect of unitarity on the pion width has already been 
discussed in [^. As it follows from our previous comments, the unitarized width Pp has a softer behaviour with p, 
which, according to our arguments in section [111 PI means a larger contribution to the conductivity as T increases, as 
compared with the non-unitarized one. 

We replace then the unitarized lAM partial waves in 0 in order to obtain the unitarized width ® and conductivity 
m- Several comments are in order here. First, we are considering only the contribution of the partial waves with 
lower angular momentum, i.e., IJ = 00,11, 20. Partial waves with J > 1 are negligible for y/s < 1 GeV 0 and so 
they are for the temperatures involved here. Second, we are not including the finite temperature corrections to the 
partial waves calculated in |^. The reason is that those are dilute gas corrections to the partial waves and therefore 
would contribute at 0{n\) to the thermal width in i.e, at the same order that other contributions that we have 
neglected. However, we should bear in mind that thermal corrections to the amplitudes describe correctly important 
features such as the changes in the mass and width of the thermal resonances, as it has been indeed observed in 
Relativistic Heavy Ion Collisions, so that their effect could be qualitatively important. Finally, we are using the same 
set of chiral parameters for the d = 4 lagrangian as in Q, ensuring that the mass and width of the p agree with the 
T = 0 values in 

The results of the unitarized conductivity are showed in Figure |S1 where we compare with the non-unitarized one 
including 0{p‘^) terms in the partial waves, i.e., t = The most relevant feature is that for temperatures 

about the pion mass, the curve for the conductivity changes notably its behaviour with respect to the perturbative 
one for T> 100 MeV. As expected from our previous comments, the unitarized conductivity decreases much more 
slowly with temperature. In fact, extrapolating this result for high enough T, the conductivity starts increasing with 
T, although very slowly as it can be seen in Figure 0 Therefore, unitarity, and ultimately the presence of resonances 
in the medium, is responsible for the change of behaviour from the low T to the high T phases. In this sense, we notice 
that the main effect arising from unitarization comes from the 00 channel which has the vacuum quantum numbers 
and is therefore more sensible to changes near the chiral phase transition. We also remark that this softe ning of the 
temperature behaviour when unitarization effects are taken into account holds also for the shear viscosity jig. 

Although clearly we should not trust our approximations much beyond T ~ 150 MeV, it is interesting to note that 
all the high-T analysis in the QGP phase shows a typical linear increasing behaviour ct oc T j^, which follows simply 
from high-T dimensional analysis. Therefore, when we include unitarity (which may not be the only relevant effect 
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T (MeV) 


FIG. 8: We show the result obtained in the dilute gas regime with the unitarized partial waves (solid line) compared to the 
perturbative ones up to 0{p^) (dashed line). 


at these temperatures), the conductivity slowly tends to its expected high-T behaviour. In fact, extrapolating our 
results to arbitrarily high temperature yields a linear behaviour for T ~ 600 MeV. 

The high-T QCD analysis are not applicable to temperatures close to the transition, say T ^ 200 MeV, where 
one should in principle appeal either to lattice results or to phenomenological models. The lattice calculation in m 
predicts a ~ 0.35T for T ~ l.STc. However, one should be careful when interpreting results coming from the lattice 
in the w ^ O’*" limit of the current-current correlators [^ . In fact, there are notable discrepancies between different 
lattice analysis. For instance, the conductivity extrapolated from the dilepton rate in for T ~ 1.5Tc,3Tc vanishes 
identically, as pointed out in |2l| . This is equivalent to state that the photon spectrum vanishes at vanishing energy 
from the timelike region (see our comments in section FV All . 

Our results show the behaviour from below the chiral phase transition, although with our present level of approxi¬ 
mation we can only make qualitative statements at this point. On the one hand, we predict a nonzero conductivity, 
which, as we discuss in the next section, is compatible with having a sizable contribution to the photon emissivity from 
the hadronic phase. On the other hand, although unitarity makes the conductivity grow with T as we approach the 
chiral phase transition (note that the nonunitarized result vanishes asymptotically) the values we obtain are still much 
smaller than those in For instance, we have a/T ~ 0.007 around T = 200 MeV. Apart from possible ambiguities 
in the lattice analysis, we should bear in mind that several effects that we have not accounted for in our perturbative 
treatment, like resummation of ladder diagrams or the contribution of heavier particles, become important at those 
temperatures. 


A. Photon production 

Another possibility to extract physical information related to the electrical conductivity is to consider the photon 
emissivity coming from the hadronic phase. The photon differential rate emerging from an equilibrated system is 
directly related with the EM current-current correlator as 

= 1 ^ 1 ) ( 30 ) 

The Ward Identity = 0 implies Po(a; 0, (f = 0) = 0 so that the electrical conductivity provides direct 

information about the vanishing energy limit (tu —> 0'*') in the static region mill simply as ujdR-y/dPq{uj ^ 0+, g = 
0) = 3Tcr(T)/(47r^), i.e, a nonzero conductivity implies a constant value for the photon yield near the origin^^- 

With this motivation in mind, we will compare our results for the conductivity for temperatures physically relevant 
in a Relativistic Heavy Ion environment with some recent theoretical and experimental analysis of the low-energy 
hadronic photon spectrum. Among the theoretical works, one can find the virial expansion approach followed in |22|| 
with nonstrange mesons only, extended in to include baryon density. More recent analysis follow the conventional 
Kinetic Theory approach, where all possible photon-producing reactions of the type 1-I-2 —> S-l-y are evaluated, 1,2,3 
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being meson degrees of freedom 7r,p, if,... parametrized in an effective lagrangian with explicit resonances included 
[r^ I 24 1^ . The present experimental data are described remarkably well for uj >200 MeV when all the relevant 
processes are evaluated and finite density effects are accounted for j25l |. However, it is worth mentioning that recent 
analysis of the WA98 experiment at CERN have produced data at lower transverse momentum, showing a slight 
systematic enhancement with respect to theoretical predictions. At those low energies, the hadronic gas contribution 
seems to dominate over the QGP one. 

The low-energy rate with nonstrange mesons in which deals then with the same degrees of freedom as our 
approach, reaches a maximum value of about 3 xl0“® GeV“^fm“® (for T =150 MeV) at w ~ 400 MeV and then 
drops to zero at the origin. When including baryons in the same approach the behaviour is the same but the 
maximum is at tu ^ 200 MeV and grows up to about 4 xl0“^ GeV“^fm“‘^. The results in for the hadronic 
gas are similar to those in p^ 1^ except that the curve with baryon density decreases monotonically for uj > 200 
MeV. In P3, more relevant channels have been included, the most interesting result being that there are also meson 
contributions that do not vanish when approaching the origin, for instance those involving -Kmi mesons. The latter 
contributions amount to almost twice more photons near the origin that in the earlier work |24| . We remark that the 
prediction in |25|| is the one used in the experimental WA98 paper to compare with data and it remains below the 
lowest energy points. Using the unitarized width, we get for the vanishing energy photon rate iTa/^ATv^) ~ 3.7 x 10“^ 
GeV^^fm”"^ at T =150 MeV. Therefore, we predict a sizable value of the hadronic photon rate at the origin. A linear 
extrapolation of the low energy curve {uj ~ 200 MeV) in to the origin gives a value of about 2 x 10“^ GeV^^fm”"^ 
at T =150 MeV. This is an indication that our prediction lies in the correct range. 

We will now try to establish a more direct comparison with experimental data, bearing in mind that our approach 
only gives information about the value very close to the origin. One then has to integrate the photon rate in CT 
over space-time and average over the photon rapidity p^ in order to obtain the experimentally measured photon 
yield ujdN^/cPq as a function of qT^ the component of the photon momentum transverse to the collision axis in 
the laboratory frame. Obviously, the results depend heavily on the hydrodynamical space-time evolution for the 
conditions applicable to a particular experiment. We will just give here rough estimates and for that purpose we will 
neglect transverse flow and assume a simple Bjorken’s hydrodynamical description PE 113 so that uj = qr cosh(j/ — 77 ) 
with y and 77 the photon and fluid rapidities respectively. Therefore, our value at a; —*■ O'*" translates directly into the 
value of the yield at qx 0+. As for the space-time integration, we use the typical value pE H^l for the nuclear 
transverse radius Ra — 1-3 fm ~ 7.7 fm for ^*^®Pb and change to proper time r and rapidity rj coordinates. In 
the Bjorken’s limit and for qx —> O’*", there is no dependence with rapidity so that we have simply to multiply by 
the expansion velocity pE ^'Hnuci = 2arccosh(-v/s/(2AGeU) ~ 10.1 for the WA98 collision energy -y/s = 158 A GeV. 
With these approximations, the photon yield at the origin is given by: 

dN^, , , , , ri 3r(r)cr(r(r)) , , , 

^-^{dT ^ -t^Ra^VtiucI j -- —Tdx (31) 

A crude estimate is obtained by assuming a purely thermal hadronic phase, i.e., a constant temperature. Taking 
proper time values Ti = 1 fm/c and r/ = 13 fm/c p3 and T = 150 MeV, this gives dN^/d^q{qx = 0) ~ 5.8 x 10^. A 
more realistic approximation is to take the cooling law of an ideal gas T{t) = Ti{TilTy^^ pE which is probably still 
rather crude for mesons at moderate or high temperatures. Taking = 3 fm/c, more appropriate for the hadronic 
phase, Tf = 13 fm/c and Ti = 170 MeV, we get Tf = 104 MeV which is of the order of the freeze-out temperature. 
Inserting this law in m with the unitarized conductivity in FiglHlwe obtain dN^/d^q{qx = 0) ~ 5.6 x 10^ although 
this result is very sensitive to variations in temperature. 

Taking the two points of smallest px in |38| and simply extrapolating them to the origin with a straight line gives 
dNj/d^q{qx = 0) ~ 5 x 10^. Therefore, our results are compatible with the recent data, in the sense of a naive linear 
extrapolation from the origin and not forgetting that we are bordering the applicability range of our approach, in 
addition to the many approximations performed to arrive to a number comparable with experiment. In any case, our 
analysis would suggest that purely thermal cuts in pion-pion annihilation (there is no T = 0 real photon production 
from TTTT annihilation for energies below 2 to^) with a nonzero pion thermal width may be a relevant effect |^. In 
this sense, it should be borne in mind that the pion thermal width is not considered in [3 EE EE Q EE since it 
does not play any role at the energies of interest for the photon yield considered in those works. On the other hand, 
our approach could be too limited to describe correctly the effect of some resonances like the oi. It is even unclear 
whether our prediction at small energies has simply to be added to previous hadronic ones, since we are extracting 
the photon rate from the self-energy and not from individual processes. What is more clear in this respect is that the 
contribution analyzed here does not come from baryonic sources. 
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VI. CONCLUSIONS 

In this work we have analyzed the DC electrical conductivity of a pion gas at low temperatures within the framework 
of Chiral Perturbation Theory. The nonperturbative nature of transport coefficients is reflected in the need of including 
the thermal pion width in the calculations, in order to avoid “pinching poles” singularities. Physically, this allows 
to account for the relevant in-medium pion collisions and gives rise to the leading order contribution in the inverse 
width consistently with Kinetic Theory. The leading order in ChPT shows a decreasing behaviour with temperature, 
behaving like cr ~ \JrriT^/T for T « TO,r as expected from the nonrelativistic limit of Kinetic Theory. For higher 
temperatures, cr ~ (m,n./T)^ for T<mT^ if unitarity corrections are not included in elastic pion scattering. 

A very important part of our analysis has been the role of higher order diagrams which, although naively suppressed, 
are enhanced by powers of the inverse width. As it happens in other theories, the dominant diagrams are uncrossed 
ladder ones, which in our case can be interpreted in terms of pion scattering of the internal lines. A careful evaluation 
of the ladder diagrams shows that they can be still considered perturbative in ChPT at low temperatures. This is 
particularly important at very low temperatures T << where the nonperturbative contributions are larger. In 
that regime we have been able to show exactly that ladder diagrams are perturbative, so that they merely renormalize 
the coefficient of the IjVT behaviour of the conductivity by subleading corrections in ChPT. Although collinearly 
enhanced, our numerical analysis shows that these corrections are typically around 5%. As temperature increases, 
ladder diagrams become more important, the most relevant contributions coming from ladders with derivative vertices. 
At temperatures near the chiral phase transition it is not clear that ladder diagrams can be neglected, since ChPT is 
not applicable. In fact, at those high temperatures we expect other effects to become important, like the presence of 
kaon states. 

Another important point concerns unitarity. We have showed that the conductivity changes qualitatively its be¬ 
haviour with temperature for T >100 MeV as a consequence of implementing unitarity of the partial waves in the 
thermal width. This provides a more physical picture as far as the behaviour of partial waves with energy and the 
presence of resonances in the thermal bath are concerned. The unitarized conductivity increases slowly for increasing 
T, which seems to be consistent with lattice and analytical analysis far beyond the transition point. In fact, the 
most important contribution in this respect comes from the partial wave corresponding to the /o(600) with vacuum 
quantum numbers / = J = 0. We have also considered some possible consequences of our results regarding the photon 
spectrum. The lack of precise experimental and lattice knowledge about the zero frequency limit, together with the 
own limitations of our approach do not allow to draw very quantitative conclusions. However, our analysis implies 
that there should be sizable effects for very low energy hadronic photon production and this result is consistent with 
recent theoretical low-energy analysis and compatible with naive extrapolations of experimental data. In our opinion, 
this makes the analysis presented here usefnl in order to understand some of the current discrepancies regarding the 
low-energy results. However, this is just but the first step and we will pursue a more detailed analysis in the near 
future, including some effects discarded here that could be important in more realistic situations, like the presence of 
kaon states, nonzero pion and baryon chemical potentials or thermal modifications of resonances. 

Summarizing, despite its limitations as the temperature approaches the phase transition, we believe that our present 
work provides a useful and physically relevant example where the diagrammatic approach to transport coefficients can 
be implemented in a controlled way. More realistic applications to physical observables such as the photon spectrum 
are under way and will be reported elsewhere. 
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APPENDIX: THERMAL CORRELATORS 


1. General Properties 

We review here some of the results needed throughout the text, mostly to fix our conventions and notation. For 
a given, rotationally invariant, correlator A(z, jgl) analytical for any complex z off the real axis, it is convenient to 
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introduce its spectral function as: 


A(z,|g1) 



dw' p{w', |(^) 
27r oj' — z 


The advanced and retarded correlators are then given by |42 


(A.l) 


l^) = ipiA(z = uj±ie, |^) 


(A.2) 


where a> € K. Note that, with our definitions p = 2ReA^ = 2Im(iA^). 

On the other hand, the imaginary-time correlator is given by the values on the discrete Matsubara frequencies: 


A,3(*a;„, l^) = A(z = Igl) (A.3) 

When performing a perturbative calculation in the imaginary-time formalism, one can express the final result as 
an analytic function off the real axis, once all the sums over internal frequencies have been performed. Then, the 
imaginary-time result is continued analytically to external real frequencies simply from (IA.2II and as: 

l^) = TiApiiuJn ^ uj±ie, |^) (A.4) 

The usual case is that when the correlator is a two-point function, like the held (j) itself or a conserved current. In 
that case one has in position space: 

= ±0(±t) < [0(<, a?), 0(0)] > 

Gf3iT,x) = < Tc0(-iT,f)0(O) > 

p(t,x) = < [0(t,f),0(O)] > (A.5) 

where the brackets denote thermal expectation values, t G R and Tq denotes time ordering along the imaginary-time 
contour C dehned as the straight line t G [0,/3]. Note that, from the above equations, p{co, |^) = —p{—uj, |g|), as 
expected for scalar spectral functions. 


2. Loop integrals 

In the analysis of the conductivity, we hnd loops involving the imaginary-time correlators: 


liitUm,\Q\]k) = iT 


tfp 




2k 


n— — (x> ' 


( 2 ^)" [ujI + W + + \p + Q? + 


(A.6) 


with fc = 0,1. 

Performing the Matsubara sum in the usual way one arrives to a function analytical in iuim off the real axis, 
so that we hnd for the retarded correlator: 

Im/^(u;, |g|;fc) = TT ^ —[I + n [siEp) + n {s 2 Ep+Q)] 

(2^) ^EpEp+Q 

x5[uj - siEp- s2Ep+q\ (A.7) 

where E'^ = \q\^ + w? and nsix) = — l) ^ is the Bose-Einstein distribution function. The two Si = S 2 

contributions require s > 4m^ (s = — |Qp) and give the unitarity cut appearing in the analysis of thermal 

scattering [^. Performing the angular integrations with the delta function, it can be written as: 
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IinJ^(c. 


\Q\\k) 


sgn(aj) 6 >(s — Aw?) 
IQtt\Q\ 


ry+ 


dy {y'^ - • 




riBiy) +nB{\uj\ -y)] 


(A. 8 ) 


where y± = (|w| ± \Q\a2p{s))/2 > m for s > Aw? and (J 2 p{s) = sjA — Aw? js is the phase space of two particles of 
equal mass m. 

The case k = 0 can be solved in terms of elementary functions: 


Im/^(a;, |( 3 |; 0 ) 


sgn(u;) 0 (s — 4m^ 
IGtt 


a’2p(s) 


2T 


■log 


1 - exp(-y+/r) 


_1 - exp(-y_/r) 


On the other hand, in the limit \Q\ O'*' (center of mass limit in the scattering case) we find: 


(A.9) 


Im/"(a;,|Q| ^0+;fc) 


sgn(w) 0 (s — 4m^) 
1071 


(Tt(w) 



(A.IO) 


where (Tt(w) = cr 2 p(oj‘^) [1 + 2nB{\oj\/2)] is the two^article thermal phase space, which can be interpreted in terms of 
enhancement and absorption in the thermal bath [ 3 - 

Let us analyze now the two si = —S 2 contributions. First, it is not difficult to see that these two integrals are 
identical and contribute only for s < 0. Unlike the unitarity cut, this cut is purely thermal, i.e, is only present for 
T 7 ^ 0. Integrating the delta function, their contribution to Im/^ is: 


Im/"(a;,|Q|;fc) 


sgn(u;) 6 »(-s) 


ther 


dy [riBiy + |w|) - riBiy)] 


IQ'k\Q\ 

X [(y2_^2)fc + ((y+ 1^1)2 


(A.ll) 


where yo = —y- = (—|i^| + |Q|o' 2 p(s ))/2 > w for s 7 ^ 0. For fc = 0: 


Im/^(w, |Q|; 0 ) 


T 


- ther 


SttIQI 


sgn(a;) 0 (—s) log 


1 _ g-lyo + li^D/T 


1 — e“2/o/'r 


(A.12) 


Now, taking the limit \Q\ ^0+ with w and T fixed implies necessarily w —> 0+ since this is a space-like contribution. 
More precisely, j/o — so that: 


Im/^(a;^0+,|Q| ^0+;/c) 

= -^0(-s)^ 

Fk 

c^ViQY;T 


- 

ther Sir IQI 



\ m T J 


where: 


(A.13) 


/ OO 

dy(y2 - rr?)^n'B{y) (A.14) 

for 0 < a; < 1. Note that Fk{x) < 0, lim^,^!- Fk{x) = 0 and that the small energy and momentum limits (I A. 1311 do 
not coincide necessarily with taking directly IQj = 0 with w 7 ^ 0 , which gives a vanishing contribution due to the step 
function. 

Following similar steps as above, it is not difficult to obtain the imaginary part of loop integrals with more powers 
of iLOn in the numerator: 


d?p 


Ja{iU}mAQ\) =iT ^ 

u=-ood 


{iuJnY 


[u)1 + \'pY + m^] {ujn + F \'P F QY + 


(A.15) 
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with a a positive integer or zero. We have analyzed the a = 0 case previously. For arbitrary a we get: 


Im (w, IQI) 

Im |Q|) 


[sgn(a;)]“ 0(s — 4m^) 

IGttIQI Jy_ 

[sgn(o;)]“+'0(-s) r 


dy{-yT [1 + riB{y) + - y)] 


(A.16) 


- ther 




'Vo 


dy[y^ + {-iny + \u:\r] 
X [nB(y + |w|) - nsiy)] 


(A.17) 
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